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1.  INTRODUCTION 


The  Geosciences  Graduate  Fellowship  program  was  established  by  the  Army  Center  for 
Geosciences  at  Coltnado  State  University  in  September,  1986.  This  five  year  program  was  funded 
by  the  Army  Research  Office  (ARO)  in  conjunction  with  the  Center  fn:  Geosciences  research 
program.  Tbe  Geosciences  Fellowship  program  served  to  increase  the  science  base  of  the  Am^ 
and  the  Nation  by  educating  students  in  technical  disciplines  while  fostering  an  q)preciati(Mi  of  the 
Army's  research  objectives,  mediods,  ami  facilities.  Tlte  coordinaticMi  between  Ae  Fellowship  and 
research  activities  provided  the  Fellows  with  an  environment  of  excellence  in  basic, 
interdisciplinary  research.  In  return,  the  Fellows  ctxitributed  dieir  talents  to  help  the  Center  meet 
its  research  go^s. 

The  Geosciences  Fellows  were  selected  on  the  basis  of  their  qualifications  demonstrated  by 
prior  academic  gradepoint.  Graduate  Record  Examinatitm  scores,  written  recommendations,  and 
their  interest  in  research  tqrics  being  perfonxMd  at  the  Cento*.  The  Fellows  were  enrolled  in 
academic  programs  while  pursuing  their  research,  leading  to  Master  of  Science  or  Doctor  of 
Philosophy  degrees.  These  degrees  were  obtained  in: 

•  Atmospheric  Science, 

•  Qvil  Engineering, 

•  Electrical  Engineering, 

•  Physics, 

•  Computer  Science, 

•  Earth  Resources,  or 

•  Psychology. 

A  total  of  18  graduate  degrees  have  been  granted  to  date  as  a  direct  result  of  the  Fellowship 
program.  Another  eight  students  still  are  woridng  to  fulfill  their  degree  requirements  following  the 
end  of  the  program  in  Novemba:,  1991.  Seven  of  the  degrees  granted  were  Ph.D.’s  and  1 1  were 
M.S.'s.  A  list  of  Fellows  and  summary  of  degrees  obtained  is  shown  in  Table  1.  We  note  that 
several  graduate  students  were  supported  by  £e  Center  for  Geosciences  research  program  as 
Graduate  Research  Assistants.  The  wmk  of  diose  students  is  not  summarized  here,  ^t  will  be 
included  in  the  Geosciences  Research  Final  Report  The  Geosciences  Fellowship  program  was 
administered  under  the  hierarchy  shown  in  Hgure  1,  whoreas  the  individual  Fellows  were 
supervised  by  their  academic  departments  acting  thrragh  a  faculty  advistv. 

This  report  contains  a  list  of  the  thesis  and  dissertaticm  titles  completed  by  the  Geosciences 
Fellows,  shown  in  Table  2.  It  also  includes  the  abstracts  of  those  treatises. 
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Table  1 


Degree  Information  for  ARO  Geosciences  Fellows 

Academic 

Fellow  Name _  Degree  Department  _ Advisor 


Raul  Alvarez 

Ph.D, 

Physics 

Dr.  J.  She 

Deborah  Anthony 

M.S. 

Earth  Resources 

Dr.  S.  Schumm 

Bonnie  Ashbum 

M.S.+ 

Atnoospheric  Science 

Dr.  E.  Reiter 

Charlotte  Atwater 

M.S. 

Atmospheric  Science 

Dr.  T.  Vonder  Haar 

James  Bossert 

Ph.D. 

Atmospheric  Science 

Drs.  W.  Cotton/R.  Pielke 

Thomas  Burke 

M.S.* 

Civil  Engineering 

Dr.  P.  Julien 

Lx>ren  Caldwell 

Ph.D.* 

Physics 

Dr.  C.  She 

Thaddeus  Cline 

M.S. 

Civil  Engineering 

Drs.  D.  Simons/P.  Julien 

Jeff  Copeland 

Ph.D.* 

Atmospheric  Science 

Drs.  W.  Cotton/R.  Pielke 

Keely  Costigan 

Ph.D.* 

Atnoospheric  Science 

Dr.  W.  Cotton 

Jeff  Fredericks 

Ph.D.+ 

Civil  Engineering 

Dr.  D.  Simons/P.  Julien 

Scot  Greenidge 

M.S.+ 

Electrical  Engineering 

Dr.  T.  Brubaker 

Benjamin  Hayes 

Ph.D.* 

Earth  Resources 

Dr.  S.  Schumm 

Russell  Huonder 

M.S. 

CcMxq)uter  Science 

Dr.  T.  Brubaker 

Andrew  Jones 

M.S. 

Ph.D.* 

Atmospheric  Science 

Dr.  T.  Vonder  Haar 

Scot  Makinen 

Ph.D.+ 

Electrical  Engineering 

Dr.  T.  Brubaker 

Robert  McCoy 

M.S. 

Physics 

Dr.  G.  Stephens 

Michael  Meyers 

M.S. 

Ph.D.* 

Atmospheric  Science 

Dr.  W.  Cotton 

Glenn  Moglen 

M.S. 

Qvil  Engineering 

Drs.  D.  Simons/P.  Julien 

Fred  Ogden 

Ph.D. 

Civil  Engineering 

Dr.  P.  Julien 

*  In  progress 

+  Left  program  without  degree 
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Table  1  (cont'd) 


Fellow  Name _  Degree 


Peter  Olsson 

M.S. 

Jerry  Richardson 

Ph.D. 

William  Smith 

M.S. 

Otto  Stein 

Ph.D. 

Joseph  Turk 

Ph.D. 

Paul  Weiler 

M.S. 

Douglas  Wesley 

Ph.D. 

Paul  Wolyn 

Ph.D.* 

Academic 

Department  _ Advisor 


Atmospheric  Science 

Dr.  S.  Cox 

Gvil  Engineering 

Dr.  P.  Julien 

Atmospheric  Science 

Dr.  S.  Cox 

Civil  Engineering 

Dr.  P.  Julien 

Electrical  Engineeting 

Dr.  V.  Bringi 

Psychology 

Dr.  R.  Loomis 

Atnx>spheric  Science 

Dr.  R.  Pielke 

Atmospheric  Science 

Dr.  T.  McKee 

*  In  progress 

+  Left  program  without  degree 
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Table  2 


Geosciences  Fellows  Thesis  and  Dissertation  Titles 


Name 

Title 

Degree 

Raul  John  Alvarez,  n 

Measurement  ttf  Trc^)Ospheric  Tenq)erature 
and  Aerosol  Extinction  Using  High  Spectral 
Resolution  Udar 

Ph.D. 

Deborah  J.  Anthony 

Stage  Dependent  Channel  Adjustments  in  a 
Meandering  River,  Fall  River,  Colorado 

M.S. 

Charlotte  A.  Atwater 

Rader-Derived  Estimates  of  Latent  Heating 

Rates  in  a  Convective  Storm 

M.S. 

James  E.  Bossert 

Regional-Scale  Flows  in  Complex  Terrain: 

An  Observational  and  Numerical  Investigation 

Ph.D. 

Thaddeus  Qine 

Development  of  a  Watershed  Information 

System  for  HEC-1  with  Application  to 

Macks  Creek,  Idaho 

M.S. 

Benjamin  Hayes 

Climatic  and  Geomorphic  Ccmttols  on  the 
Suspended  Sediment  and  Salt  Loads  in  the 

Upper  Colorado  River  Basin 

Ph.D. 

Russell  J.  Huonder 

Visualization  of  Three-Dimensional  Scientific 

Data  Using  Octrees 

M.S. 

Andrew  S.  Jones 

Nficrowave  Remote  Sensing  of  Qoud  Liquid 

Water  and  Surface  Emittance  Over  Land  Regitxis 

M.S. 

Robert  F,  McCoy,  Jr. 

Frequency  Chirp  of  a  Hybrid  Doppler  Lidar 

M.S. 

Michael  P.  Meyers 

An  Evaluation  of  the  Factors  Affecting  Wintertime 
Quantitative  Precipiuition  Forecasts  in  an 

Explicit  Qoud  Model  over  Mountainous  Terrain 

M.S. 

Glenn  E.  Moglen 

The  Effects  of  Spatial  Variability  of  Overiand 

Flow  Parameters  on  Runoff  Hydrogrtqphs 

M.S. 

Fred  L.  Ogden 

Two-Dimensional  Runoff  Modeling  with 

Weather  Radar  Data 

Ph.D. 

Peter  Olsson 

Inference  of  Horizontal  Teiiqterature  Gradients 
using  Passive  Radiometric  MetiKxls 

M.S. 

Jerry  R.  Richardson 

The  Effect  of  Moving  Rainstnms  on  Overiand 

Flow  Using  One-Dimensional  Finite  Elements 

Ph.D. 
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Name 

Table  2  (cont'd) 

Title 

Degree 

William  L.  Smith,  Jr. 

The  Broadband  Radiative  Properties  of  Cirrus 
Qouds  Deduced  from  Aircraft  Measurements 
During  FIRE 

M.S. 

Otto  R.  Stein 

Mechanics  of  Headcut  Migration  in  Rills 

Ph.D. 

Joseph  Turk 

Microwave  Remote  Sensing  of  a  Precipitating 
Atnx)sphere 

Ph.D. 

Paul  A.  Weiler 

Social  Loafing  and  Vigilance  Decrement  While 
Monitoring  E^erent  Levels  of  Simulated  Air 
Traffic 

Ph.D. 

Douglas  A.  Wesley 

The  Influence  of  Topography  on  Colorado 

Front  Range  Snowstorms 

Ph.D. 

5 


Hgure  1.  Administration  of  Geosciences  Fellowsh^  Program 
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2 .  ABSTRACTS  OF  THESES  AND  DISSERTATIONS 


ABSTRACT 

MEASUREMENT  OF  TROPOSPHERIC  TEMPERATURE  AND  AEROSOL 
EXTINCTION  HIGH  SPECTRAL  RESOLUTION  LIDAR 

Raul  John  Alvarez,  n 
Physics  Department 


A  high  spectral  resolution  Rayleigh-Mie  lidar  system  capaUe  of  measuring  relative  ta]:q)erature 
and  aerosol  profiles  in  the  troposphere  has  beat  dt^elopcd.  The  use  of  a  narrowband  laser 
transmitter  sdong  with  an  atomic  vt^xn*  cell  as  a  variable  width  bandstop  filter  allows  for  the 
separation  of  molecular  and  aerosol  scattering  signals.  Using  r.vo  of  these  filters  with  different 
widths  and  a  newly  develc^ied  lidar  inversion  mediod  it  is  possible  to  detennine  relative 
tetiq)mture,  pressure,  and  density  profiles  of  dte  tr(^x>sphere.  This  separation  also  makes  a 
quantitative  measurement  of  aerosol  extinction  coefficient  profiles  possible. 

In  principle,  the  information  required  in  addition  to  the  lidar  signals  fix*  determining 
atmospheric  and  aerostd  prop^es  with  this  Rayleigh-Mie  lidar  is:  the  theoretical  Rayleigh- 
Brillouin  spectrum  for  scattering  of  narrow  bandwiddi  laser  light,  independent  measurements  of 
the  filter  transmission  functions,  the  assun^tion  that  the  atmos^cre  acts  as  an  ideal  gas  in 
hydrostatic  equilibrium,  and  an  approximate  atmospheric  pressure  at  one  reference  altitude. 

Profiles  of  the  tropospl^c  state  variables  tenq)mture,  pressure,  and  density  are  determined 

from  signals  averaged  over  20  minutes.  Once  profiles  of  die  atmospheric  variables  are  determined, 
the  aerosol  extinction  coefficient  profiles  are  oSknilated. 

A  systematic  offset  in  the  lidar  tenqterature  profites  is  seen  however  in  practice  when  conqiared 
to  balloon  sonde  measurements.  Further  labcxatray  experiments  have  established  that  the 
systematic  offset  is  real  and  is  consistent  between  labcmtory  and  field  measurements.  This 
repeatable  offset  can  be  eliminated  by  normalizing  the  tenq)erature  at  a  reference  altitude  to  an 
independently  measured  or  estimated  value. 

The  development  and  in:q)lementation  oi  this  high  spectral  resolution  lidar  and  the  first 
simultaneous  measurements  of  trcqx)spheric  temperature  and  aerosol  profiles  from  0.2  km  to  8  km 
are  presented.  The  precision  at  1  Im  dtitude  is  ±10  K  for  temperature  measurement,  ±3%  for 
backscatter  ratio,  and  ±5%  fcx  aerosol  extinction  coefficient  At  5  km,  these  values  are  ±12  K  for 
temperature,  ±4%  f<x  backscatter  ratio,  and  ±100%  for  extinction  coefficient  Suggestions  to 
determine  the  cause  of  the  offset  and  to  reduce  the  relative  uncertainties  are  also  discussed. 
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ABSTRACT 


STAGE  DEPENDENT  CHANNEL  ADJUSTMENTS  IN  A 
MEANDERING  RIVER,  FALL  RIVER,  COLORADO 

Deborah  J.  Anthony 
Earth  Resources  Depaiunent 


Dam  faUure  increased  bedload  supply  to  Fall  River,  a  very  sinuous,  low  gradient,  snowmelt- 
fed  stiea^  by  a  factor  of  1000.  This  setfiment  is  mobile  at  nt^uly  all  discharges,  and  thus  channel 
topography  is  able  to  change  with  stage.  In  a  two-bend  study  reach,  the  relationships  between 
tl<w  and  mobile  bed  material  were  studied  by  measuring  bedoad  transport  rates,  bed  and  water 
surface  topography,  and  downstream  and  cross-stream  velocity  patterns  on  22  cross  sections 
Measurements  were  made  at  different  discharge  levels  through  the  1986  runoff  season. 

At  bankfull  discharge,  bend  cross  sections  showed  maximum  asymmetry,  witfi  point  bar 
platfcnms  built  up  to  the  water  surface  and  thalwegs  excavated  to  their  greatest  depth,  ^ss- 
stream  flows  included  zones  of  high  turbulence,  containing  bed  and  bank  cells  in  the  thalweg 
region,  with  helical  flow  occurring  over  the  point  bar  slope  and  outward  flow  only  occurring  over 
the  pomt  bar  platforms.  B^oad  transport  was  greatest  over  the  point  bar  slope,  and  almost  zero 
through  the  thalweg.  Sorting  was  continuous  drrough  each  bend. 

At  intCTm^ate  flows,  cross  section  symmetry  increased  due  to  lateral  erositm  of  the  point  bar 
platform.  Helical  flows  occupied  a  larger  portion  of  individual  cross  sections,  inrlnriing  the 
thalweg.  In  planform,  however,  the  zone  of  helical  flow  had  ^Ut  into  two  smaller  cells  in  each 
bend.  Be^oad  u^sport  had  shifted  outward  from  its  high  flow  position,  due  to  point  bar 
ptogradation  towards  the  cutbank,  so  that  some  sediment  movement  was  recorded  in  the  thalweg. 

A  flows,  characteristic  of  fall  and  winter,  cross  secti(»  symmetry  was  at  a  maTimpn^ 

Actional  erosion  of  the  high  point  bars  had  occurred,  but  the  greatest  change  was  due  to  thalweg 
filhng.  Cross-stream  flows  were  weak  and  chaotic  for  the  most  part,  and  bedload  transport  was 
more  uniform  across  the  channel  than  at  higher  discharges.  Sorting  was  discontinuous,  and  it  was 
confined  to  small  pools  in  the  otherwise  flat  channel 

In  this  mobile  bed  channel,  bed  topography  is  a  function  of  cross-stream  flow  and  bedload 
transport  patterns,  which  are  control!^  by  stage. 


8 


ABSTRACT 


RADAR-DERIVED  ESTIMATES  OF  LATENT  HEATING  RATES 
IN  A  CONVECTIVE  STORM 

Charlotte  A.  Atwater 
Atmospheric  Science  Department 


The  importance  of  latent  energy  in  sttMm  d'mamics  and  the  global  energy  balance  has  long  been 
recognized  by  atmospheric  researchers.  This  study  focuses  on  the  latent  heating  and  cooling  in  an 
isolated,  microburst-producing  thunderstorm  in  Alabama  on  the  20th  of  July,  1986.  Data  used  in 
this  work  were  record  by  NCAR  CP-2  radar.  10  cm  horizontal  and  vertic^  reflectivity  signals 
(Z//  and  Zv,  respectively)  are  used  to  calculate  Ae  difference  reflectivity,  Z^p  =  101og(Z//  -  Zy), 
which  is  Aen  used  to  ctdculate  Ae  fraction  of  Ae  horizontal  reflectivity  signal  arismg  fiW  ice 
hydrometeors.  This  information  is  Aen  used  to  quantify  liquid  and  ice  water  contents  and  Ae 
rainfall  rate.  The  time  rate  of  change  of  Aese  thrw  quantities  arc  then  used  to  cterive  latent  heating 
rates. 

The  fiction  of  Ae  reflectivity  signal  arising  frc»n  ice  hydrometeors  is  based  on  Ae  degree  of 
decorrelation  between  Ae  Zqp  and  ^  signals  m  mixed-phase  regions.  Liquid  water  contents  arc 
calculated  usmg  a  combination  of  1)  an  empirical  relation  relating  liquid  water  content  to  radar 
reflectivity  for  mixed-phase  regions  and  2)  calculating  liquid  water  ccmtents  from  radar-derived 
parameters  of  Ae  raind^  size  Astribution  fw  pure  liquid  water  regions.  Ice  water  contents  are 
calculated  from  an  empirical  relation  relating  liquid  water  content  to  radar  reflectivity  wi  A 
conations  made  to  account  for  differences  between  Ae  iiKlices  of  rcfintetion  and  typical  densities 
of  liquid  and  ice  hydrometeors.  Finally,  rainfall  rates  are  calculated  usmg  an  empirical  relation 
relating  rainfall  rate  to  coincident  differential  reflectivity,  Zop,  and  Z  signals.  Estimated  enor 
levels,  including  mstrument  error,  arc  approximately  ±10%  for  Ae  time  rate  of  change  of  liquid 
water  content,  ±70%  for  Ae  time  rate  of  change  of  ice  water  content,  and  ±40%  for  rainfall  rate. 

Calculated  vertical  ice  fraction  profiles  indicate  that  Ae  stram  ctmsisted  almost  entirely  of  liquid 
water  during  its  early  stages.  Glaciation  began  at  Ae  top  of  Ae  storm  and  continued  to  progress 
downward  through  Ae  stonn  wiA  Ae  descent  of  Ae  core  of  maximum  reflectivity  signd.  "nie  ice 
fraction  profile  increased  sli^tly  from  nearly  0.0  to  about  0. 1  at  Ae  surface  at  the  time  the 
microburst  was  observed  to  impact  Ae  ground.  This  increase  is  most  likely  a  result  of  pea-sized 
hail  reachmg  Ae  surface.  These  conclusions  confirm  Aose  drawn  in  oAer  stuAes  of  Ae  20  JAy 
1986  case. 

The  latent  heating  resAts  mAcate  that  heating  and  cooling  caused  by  condensation  and 
evaporation,  rcs^)ectively,  account  for  tougAy  90%  of  Ae  totA  heating  and  cooling  m  this  storm, 
with  the  rcmaiiting  10%  cotmng  from  melting  and  freezing  processes.  Net  latent  toting  rates 
(totA  heating  rate  minus  totA  coolmg  rate)  mcreased  nearly  steadily  from  Ae  early  stages  of  Ae 
storm  unA  about  10  minutes  before  Ae  rmcroburst  occurred,  at  which  point  it  began  to  decrease. 
The  net  latent  heating  rate  went  from  positive  to  negative  less  Aan  5  minutes  before  Ae  rmcroburst 
impacted  Ae  ground.  Towards  Ae  end  of  Ae  decaying  stage  of  Ae  stonn,  net  latent  heating  was 
approaching  zero. 
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ABSTRACT 


REGIONAL-SCALE  FLOWS  IN  COMPLEX  TERRAIN: 

AN  OBSERVATIONAL  AND  NUMERICAL  INVESTIGATION 


James  E.  Bossert 
Atmospheric  Scimce  Dq)artment 


An  observational  program  has  been  conducted  to  obtain  information  concerning  thermally- 
driven  flows  in  complex  terrain  on  meso-^  to  meso-a  scales  (100  -  500  km).  Data  were  collected 
from  remote  surface  observing  systems  at  exposed  mountaintop  locations  throughout  the  state  of 
Colorado,  over  the  summers  of  1984-1988.  These  field  e:q)eriments  have  been  called  the  Rocky 
Mountain  Peaks  Experiments  (ROMPEX).  The  observations  from  ROMPEX  have  been 
supplemented  with  data  from  other  remote  surface  networks,  special  soundings,  upper-air 
observations,  and  radar  and  lightning  strike  information  to  provide  an  adequate  description  of  the 
flows  and  weather  of  interest 

The  observations  have  shown  the  development  of  a  recurrent  ’’regional-scale”  circulation 
system  across  the  Colorado  mountain  barrier,  (^)erating  on  a  diurnal  time  scale.  The  basic 
structure  of  the  flow  system  consists  of  a  daytime  inflow  phase  toward  the  mountains  altxig  the 
Continental  Divide,  and  a  nocturnal  outflow  away  from  diis  high  terrain.  Long-term  averages 
show  this  circulation  system  to  be  the  dominant  wiiKl  pattom  at  several  high  altitude  stations, 
reve^g  its  climatological  significance.  Attention  has  been  focused  upon  the  nocturnal  phase  of 
the  circidation  system  along  die  western  slope  of  the  nx>untain  barrier.  Here,  the  winds  are 
particularly  strong  and  from  a  southeasterly  direction,  which  is  generally  counter  to  die  upper-level 
winds,  and  onset  abrupdy  in  early  evening  with  steady  flow  thereafter.  Soundings  have  shown 
this  nocturnal  current  to  be  shallow  and  within  a  distinct  stable  air  mass.  Clonvective  stmas  are 
found  to  enhance  this  southeasterly  flow  regime. 

Numerical  simulations  have  been  performed  widi  the  Colorado  State  University  Regitxial 
Atmospheric  Modelling  System  (CSU-RAMS)  to  provide  further  insight  into  die  physical 
mect^is^  forcing  the  observed  re^onal-scale  cuculation  systeiiL  ’^e  model  simulations  include 
both  idealized  two-  and  three-dimensional  experiments,  as  well  as  a  three-dimensional  case  study 
experiixMnt  using  actual  data  fix'  the  initialization.  The  three-dimensional  simulations  use  two-way 
interactive  grid  nesting  and  a  realistic  representation  of  tqiography  over  the  region  of  interest 

The  idealized  three-dimensional  experiment  showed  that  thermal  forcing  over  realistic 
topography  in  conditions  of  negligible,  ot  weak  ambient  flow,  is  capable  of  producing  many  of  the 
flow  features  observed  throughout  the  diurnal  cycle.  This  experiment  further  showed  how  the 
deep  mountain-plains  solenoid  altxig  and  above  the  Born  Range  crest  evolves  in  late  aftorKxxi  into 
a  shallower  density  current,  which  tiiien  pn^agates  westward  over  the  mountains  of  the  western 
slope.  This  unexpiroted  flow  phenomenon  is  ^  primary  process  responsible  for  the  strong 
nocturnal  southeasterly  winds  found  in  observations.  Sensitivity  experiments  show  that  the 
particular  terrain  configuration  through  an  east-west  cross-section  of  tiie  Colorado  mountains  is 
important  to  the  generation  of  this  unusual  circulation.  The  strong  thomal  gradient  {xoduced  by 
diffnential  heating  of  the  topography  is  the  primary  driving  fence  in  the  densiQr  current  evolution. 
(Coriolis  influence  maintains  the  steady  nocturnal  south-southeast  winds  over  die  western  slope. 
Additional  experiments  show  that  the  diumally  evolving  regional-scale  circulation  system  over  the 
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Colorado  Rocky  Mountains  is  a  robust  feature  which  can  occur  over  a  range  of  ambient  flow  and 
stratification  conditions.  Soil  moisture  experiments  reveal  that  wet  soil  altxig  the  eastern  slq)e  and 
dry  along  the  western  slope  aids  the  development  of  the  westward  prc^agating  density  current 

The  diurnal  evolution  of  the  circulation  system  on  die  case  study  day  was  in  fair  agreement 
with  many  of  the  observed  circulation  features.  This  experiment  also  revealed  that  synoptic-scale 
forcing  can  influence  the  development  of  the  regional-scale  circulaticais  in  preferential  regitxts 
along  the  eastern  slope  of  the  mountain  barrier.  As  a  result  of  the  numeric^  experiments  four 
phases  of  the  thermally  forced  regional-scale  diurnal  circulatitHi  system  have  bmn  identified. 

These  consist  of  a  daytime  mountain  boundary  layer  development  phase,  a  late  afternoon 
transitional  phase,  an  evening  propagating  density  current  phase,  and  a  late  night  adjustment  phase. 
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ABSTRACT 


DEVELOPMENT  OF  A  WATERSHED  INFORMATION  SYSTEM  FOR  HEC-1 
WITH  APPLICATION  TO  MACKS  CREEK,  IDAHO 

Thaddeus  Cline 
Civil  Engineering  Department 


A  micro  computer  based  Watershed  Information  System  (W.I.S.)  is  developed  to  assist  in  the 
preparation  of  HEC-1  input  files.  This  system  consists  of  thiM  phases.  Phase  I  utilizes  the 
capabilities  of  AutoCAD  version  9  and  three  AutoLISP  programs:  BASINS,  PLANES  and 
CHANNELS,  to  extract,  organize  and  display  watershed  ^ta.  Phase  n  uses  a  FORTRAN 
program:  CN,  to  calculate  some  HEC-1  parameter  values.  Htase  m  utilizes  a  PASCAL  program: 
HECUPDATE,  to  create  HEC-1  input  Bles.  Input  into  the  system  includes  topographic,  soils, 
land  use,  watershed  geometry  data,  and  a  ^elet^  HEC-1  input  file.  Ou^ut  from  the  system 
consists  of  a  summary  User  Reference  File,  a  Soils  File,  a  I^d  Use  I^e,  a  watershed  Geometry 
File,  a  Curve  Number  File  and  HEC-1  input  file  which  is  ready  to  run.  TTiis  thesis  presents  the 
W.I.S.:  describing  its  use,  applying  it  and  providing  program  listings. 

The  W.I.S.  is  applied  in  a  HEC- 1  nsodeling  study  of  Macks  Creek  Watershed:  a  small,  arid 
mountainous  watershed  located  in  southwest  IdalK>.  First  the  W.I.S.  generates  a  database  and 
HEC-1  input  files.  Then  these  files  are  used  to  calibrate  the  model  to  a  rainfall-runoff  event  on 
June  5  and  6, 1967.  The  input  files  from  the  calibrated  model  are  then  used  to  simulate  the 
response  of  the  watershed  to  a  heavy  rainfall  event  on  August  23, 1965.  The  W.I.S.  develop 
and  appUed  in  this  thesis  greatly  facilitates  watershed  data  extractitm  and  organizadtxi,  data  display, 
HEC-1  parameter  calculation  and  the  creation  of  HEC-1  input  files.  Two  FORTRAN  programs 
HECQ.FOR  and  HEC123.FOR  refonnat  HEC-1  output  into  files  which  can  be  import  by  a 
spreadsheet  program  for  graphical  representation  of  simulation  results.  HECQ.FOR  also  compares 
the  simulated  and  the  measi^  hydrogrtphs  by  calculating  absolute  and  relative  error,  the  mean 
and  standard  deviation  of  these  errors  and  the  volume  of  tlto  difference  between  the  simulated  and 
the  measured  hydrographs. 

It  is  shown  that  using  abnormally  low  values  for  roughness  coefficients  in  HEC-1  had  the 
same  effect  on  the  simulated  hydrognph  as  the  introduction  of  collector  channels.  It  is  also  shown 
that  using  the  base  flow  correction  option  of  HEC-1  modifies  output  hydrographs  in  a  way  that  can 
erroneously  indicate  that  volume  of  direct  runoff  is  not  preserved  when  varying  routing 
parameters.  This  is  due  to  a  "masking"  effect  justifying  caution  in  hydrograph  analysis  using  this 
option.  Finally,  a  sensitivity  analysis  of  some  EEC-1  parameters  is  perfonned.  Results  of  &s 
analysis  are  useful  for  developing  and  operating  the  W.LS.  and  for  o^el  caUlnation.  The 
sensitivity  analysis  shows  that  HEC-1  is  very  sensitive  to  values  of  curve  number  and  initial 
abstraction,  moderately  sensitive  to  values  of  overland  flow  roughness  coefficient  and  overland 
flow  length,  and  slightiy  sensitive  to  open  channel  flow  roughness  coefficient  and  ccxnputational 
time  interval. 
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VISUALIZATION  OF  THREE-DIMENSIONAL  SCIENTIFIC  DATA 

USING  OCTREES 

Russell  J.  HuondCT 
Computer  Sciences  Department 


Meteorological  senses  with  three-dimensional  distal  output  have  created  an  interest  in 
graphical  rendering  of  atmospheric  condidcxis  Sex'  scientific  visualizaticai.  The  goal  of  visualizatitm 
is  improved  understanding  of  atmospheric  phentmaena.  The  renderings  should  be  fast,  and  should 
accurately  represent  the  data.  A  traditional  approach  to  three-dimensitmal  object  rendering  is  to 
display  the  ^ta  using  a  geometric  modeling  system,  but  because  these  systems  model  object 
surfaces,  volumetric  data  sets  can  not  be  effectively  represented.  The  snoooth  surfaces  these 
systems  produce  are  also  undesirable  since  data  integrity  can  be  lost. 

In  this  thesis,  the  octree  structure  is  presented  as  a  data  structure  which  efficiently  and 
effectively  represents  the  observational  data.  An  introduction  to  this  structure  is  presen^  alcmg 
with  advantages  and  disadvantages.  Observatitmal  data  sets  used  are  discussed  dong  with  a 
description  of  how  the  data  is  encoded.  Various  existing  algorithms  for  displaying  the  stracture  are 
reviewed,  and  then  a  new  algorithm  for  visualization  of  the  octree  using  arbitrary  rotations  is 
presented.  Within  this  discussion  various  illumination  models  are  considered,  as  well  as  other 
parameters  which  permit  a  focused  view  of  the  data  such  as  slicing  and  thresholding.  A  technique 
for  fusing  data  sets  is  then  presented  with  recommendations  for  future  work  in  this  area. 

The  algorithms  devel(^)ed  for  manipulation  of  the  octree  allow  three-dimensiond  visualizatitm 
of  scientific  data.  Arbitrary  rotation,  data  fusion,  thresholding,  slicing,  and  the  various 
illumination  models  are  demonstrated  using  infrared  satellite  dto  and  radar  data.  The  results  are 
illustrated  via  multiple  photographs  of  octrees. 
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MICROWAVE  REMOTE  SENSING  OF  CLOUD  LIQUID  WATER  AND  SURFACE 

EMITTANCE  OVER  LAND  REGIONS 

Andrew  S.  Jones 
Atmospheric  Science  Department 


Microwave  remote  sensing  of  cloud  liquid  water  has  largely  been  limited  to  areas  over  ocean 
surfaces.  This  study  uses  data  from  a  new  microwave  instrument,  the  SSM/I  on  a  polar-orbiting 
DMSP  satellite,  and  infrared  and  visible  data  from  die  VISSR  instrument  on  the  GOES  satellite  in 
geostationary  orbit  The  region  selected  for  the  study  was  an  area  of  SOO  km  x  SOO  km  centered  on 
northeast  Colorado  during  die  Erst  week  of  August  1987.  The  SSM/I  instrument  has  new  high 
frequency  channels  (85.5  GHz)  which  are  more  strongly  attenuated  by  cloud  liquid  water  than 
channels  on  previous  instruments.  This  allows  fc^  the  estimation  of  integrated  cloud  liquid  water 
based  on  the  microwave  brightness  temperature  depression  caused  by  attenuation  and  emissicm  of 
microwave  radiation  at  the  colder  cloud  levels.  Atmospheric  attenuatitMi  due  to  oxygen  and  water 
vapor  is  determined  using  a  millimeter-wave  propagation  model  (MPM).  The  Rayleigh 
approximation  is  used  for  the  calculatitxi  of  cloud  liquid  water  attenuation. 

Surface  emittance  measurements  at  the  SSM/I  frequencies  were  made  with  the  aid  of  co-located 
GOES  infrared  data  during  clear  sky  ctmditions.  Images  produced  of  the  retrieved  surface 
emittances  suggest  a  strong  influence  by  wet  surfaces  caused  by  precipitation  and  irrigation.  Error 
analysis  results  indicate  absolute  erras  of  ±0.012  for  surface  emittance  retrievals  for  the  85.5  GHz 
channels. 

Integrated  cloud  liquid  water  retrievals  show  good  qualitative  agreement  with  other  available 
data  sources.  Numerical  error  sensitivity  analysis  and  conqrarison  of  integrated  cloud  liquid  water 
retrievals  for  the  vertical  and  horizcmtal  polarizaticxis  show  error  estimates  of  0. 15  kg-m-^ 
including  instrument  noise.  A  bias  between  the  horizontal  and  vertical  polarizations  of  the  85.5 
GHz  channels  was  noticed  in  the  retrieved  integrated  cloud  liquid  water  amounts.  The  bias  appears 
to  ^  due  to  a  relative  instrument  errOT  between  channels  of  approximately  1.5  K.  Absolute  error 
estimates  of  the  integrated  cloud  liquid  water  retrievals  are  unavailable  but  calibration  of  the  method 
should  be  possible  if  quantitative  integrated  cloud  liquid  water  amounts  are  known. 
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FREQUENCY  CHIRP  OF  A  HYBRID  DOPPLER  LIDAR 

Robert  F.  McCoy  Jr. 

Physics  Department 


A  hybrid  Doppler  lidar  utilizes  a  pulsed  infrared  carbon  dioxide  laser  as  the  transmitter  in  a 
laser  radar  system  for  measuring  wind  velocity  and  direction.  Some  of  the  laser  light  is 
backscatteied  by  aerosols  that  are  carried  by  the  wind.  The  scattered  light  is  shifted  by  an  amount 
proportional  to  the  aerosols'  radial  velocity.  If  the  tiequency  of  the  transmitting  laser  changes 
throughout  a  pulse,  the  velocity  estimate  of  the  wind  field  is  affected. 

The  primary  laser  within  the  system  is  the  hybrid  cavity  composed  of  two  sections.  The  first 
section  is  a  continuous-wave  (cw)  laser,  and  the  second  is  a  transverse  excited  annospheric 
pressure  (TEA)  laser  which  is  basically  a  high-power  pulsed  amplifier.  In  the  past,  (^02  Doppler 
lidars  used  custom-built  TEA  sectitms  costing  between  one-quarter  and  one-half  million  dollars. 
The  lidar  being  built  for  the  Center  for  Geosciences  at  Colorado  State  University  utilizes 
commercial  off-the-shelf  components  which  result  in  an  order-of-magnitude  reduction  in  the  cost 
of  the  TEA  section.  The  commercial  TEA  sections  do  not  appear  to  introduce  more  problems  than 
a  custom-built  TEA;  their  problems  and  solutions  are  addressed  in  the  thesis. 

The  TEA  section  is  a  gas  laser  and  exhibits  frequency  chirping  during  a  pulse.  When  the 
density  of  the  gas  within  the  TEA  section  changes,  the  ft^uency  of  the  hybnd  cavity  radiation  can 
also  change.  To  minimize  the  uncertainty  in  the  velocity  estimate,  the  frequency  chi^  should  be 
reduced  to  a  minimal  value  and  the  resulting  chirp  must  be  known.  Three  different  techniques  to 
characterize  the  chirp  are  presented,  along  with  methods  to  reduce  die  chirp.  The  chirp  is  measured 
directly  and  parameters  that  affect  the  chi^  are  identified.  Two  of  the  experiments  yield  the  chirp 
toward  the  end  of  the  pulse  where  the  direct  measurement  techniques  tend  to  provide  poor  data. 
They  also  provide  information  on  the  optical  quality  of  the  cavity  before,  during,  and  after  the 
pulse.  Suggestions  fn*  modifying  parameters  to  c^timize  the  hybrid  cavity  design,  thus  reducing 
the  chirp,  include  increased  cavity  length,  larger  niode  volumes,  vibratirai^  isolation  of  TEA  laser 
and  shorter  pulse  lengths. 
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AN  EVALUATION  OF  THE  FACTORS  AFFECTING  WINTERTIME 
QUANTITATIVE  PRECIPITATION  FORECASTS  IN  AN  EXPLICIT  CLOUD 
MODEL  OVER  MOUNTAINOUS  TERRAIN 

Michael  P.  Meyers 
Atmospheric  Science  Department 


A  prolonged  orographic  precipitation  event,  in  a  relatively  steady-state  synoptic  envirMiment 
occurred  over  the  Sierra  Nevada  Range  in  centi^  Califonia  on  12-13  February  1986.  This  well 
documented  case  was  investigated  via  the  non-hydrostatic  version  of  the  CSU  Regional 
Atmospheric  Modeling  System  (RAMS).  Two-dimensional,  cross  barrier  simulations  have  shown 
the  feasibility  of  producing  a  quantitative  precipitatitm  fraecast  (QPF)  with  an  explicit  cloud  model. 
The  simulate  flow  fields,  microphysical  stmcture,  and  precipitation  distribution  con^ared  quite 
well  with  observations.  TTie  experiment,  however,  exhibited  a  profound  sensitivity  to  the  input 
sounding.  Initializing  with  a  sounding  which  is  representative  of  the  upstream  environment  is 
most  critical  to  the  success  of  the  simulation. 

The  flow  fields  showed  little  sensitivity  to  the  model  microphysics.  Sensitivity  to 
microphysical  parameters  have  shown  that  the  removal  of  a  precipitation  process,  such  as  graupel 
or  second^  ice  production  due  to  riming,  is  con^nsated  by  the  enhancement  of  other  processes 
such  as  rain  and  aggregation.  The  sensitivity  of  lowering  the  graupel  density  effectively  changed 
the  graupel  sf^es  to  a  graupel-aggregate  hybrid  which  dominated  the  ice-phase  precipitation 
budget,  resulting  in  a  broadening  of  tlw  precipitaticm  distribution  on  the  lower  half  of  the  barrier. 

The  flow  fields  were  sensitive  to  the  absence  of  the  Oxriolis  term,  with  the  u-component  winds 
weaker  at  the  base  of  the  barrier.  The  precipitation  distribution  was  slightly  affected  with  Coriolis 
turned  off,  with  a  mote  pronounced  nM-mountain  peak  than  in  the  control  run,  which  resulted 
from  enhanced  vertical  velocities  over  the  txuiier.  The  Coastal  Range  produced  a  seeder-type 
cloud  but  enhancement  of  precipitatitxi  was  not  evident  Sensitivity  to  the  diabatic  influences  of 
melting  were  seen  in  the  microphysical  fields  but  showed  little  influence  on  tiie  flow  fields  and 
resultant  precipitation  distribution.  A  standing  15  km  wave  was  evident  upstream  of  the  barrier 
crest  and  its  origin  was  not  sensitive  to  latent  heat  effects  or  the  Coastal  Rwge. 
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THE  EFFECTS  OF  SPATIAL  VARUBILITY  OF 
OVERLAND  FLOW  PARAMETERS  ON  RUNOFF  HYDROGRAPHS 


Glenn  E.  Moglen 
Civil  Engineering  Itepartment 


.5®  of  overland  flow  system  depends  largely  on  the  degree  of  equilibrium 

e>mibited  by  the  system.  Large  variatitMis  in  discharge  result  from  partial  equilibrium  conditions, 
while  systems  m  complete  equilibrium  produce  only  small  variations  in  discharge.  The  use  of 
spati^ly  averaged  values  to  quptify  hydrologic  parameters  for  overland  flow  noodeling  is, 
uwrefore,  msufficient  information  to  axxurately  simulate  discharge  hydrographs  when  the  system 
tong  model^  has  not  reached  equilibrium.  In  contrast,  a  system  at  equilibrium  exhibits  a  peak 
discharge  mdependent  of  the  spatial  distribution  of  the  watershed  char^teristics. 

Uncorrelated  spatial  variability  in  each  of  the  input  parameters  causes  differing  relative 
variability  in  overland  flow  dischvge.  OverlaiKl  flow  ^scharge  is  noost  sensitive  to  spatial 
v^ations  in  excess  rainfall  intensity.  Manning’s  "n"  and  width  produce  comparable  variability  in 
disch^ge.  Vmations  in  slcqie  have  the  smallest  effect  on  overland  flow  discharge.  These  results 
are  relauvely  independent  of  the  spatial  distribution  functions  used  for  the  perturbations. 

‘Hie  procedures  developed  to  analyze  effects  of  uncorrelated  spatially  varied  parameters  are 
appued  using  spatially  correlated  slope  input  data.  This  is  done  using  the  spectral  distribution  of 
terrain  height  values  for  a  mountainous  region  near  Steamboat  Springs,  Colorado.  These  spectra 
^  used  to  v^  the  slope  parameter  in  the  Manning  equatioa  The  results  corroborate  the  earlier 
tindmgs  regarding  the  denxinstrated  variability  of  overland  flow  discharge  at  various  degrees  of 
equihbn^  The  results  also  indicate  that  spatial  correlation  of  a  parameter  has  an  influence  on  tiie 
magmtude  and  timing  of  the  discharge  variations. 

TTie  results  provide  a  rneans  of  determining  an  appropriate  grid  spacing  for  sufficiently 
d^cnbing  variability  Ain  a  watershed.  For  durations  of  rainfall  greater  Aan  equilibrium,  the 
effeas  of  spatial  v^bility  of  overland  flow  parameters  ate  small.  Rearrangement  of  the  equation 
for  the  time  to  equilibrium  yields  an  expression  for  the  desired  grid  spacing.  This  grid  spacing 
l^vidcs  a  incaiu  of  ensuring  thut  complete  equihbnum  will  be  achieved  within  all  grid  areas, 
therefore  mimmizing  the  error  due  to  partial  equilibrium  variations. 


ABSTRACT 


TWO-DIMENSIONAL  RUNOFF  MODELING  WITH  WEATHER  RADAR  DATA 

Fred  L.  Ogden 
QvU  Engineering  Department 


The  objective  of  this  study  is  to  examine  the  effea  of  precipitadtm  data  spatial  and  temporal 
resolution  on  the  surface  runoff  calculated  by  one-dimensional  and  two-dimensional,  physically 
based,  distributed  parameter  runoff  models. 

The  model  sensitivity  to  rainfall  data  ten^)Ofal  resoluticm  is  examined  on  both  (xie-dimensi(xial 
and  two-dimensional  impervious  surfaces.  Studies  are  performed  using  a  Monte  Carlo  approach, 
where  fifty  simulations  are  performed  fcx'  each  value  of  tempcmd  sampl^g  resoluticm  and  rainfall 
duration.  Results  indicate  that  the  sensitivity  of  both  one-dimensional  and  two-dimensional  runoff 
geometries  increases  with  increasing  temptxal  sampling  resolution  and  with  increasing  rainfall 
duration.  The  sensitivity  reaches  an  asymptote  as  Ae  rmnfall  duration  exceeds  the  time  to 
equilibrium,  defined  by  the  average  rainfall  intensity.  The  asymptotic  value  of  relative  sensitivity 
increases  with  the  square  root  of  tiie  temporal  sampling  interval. 

The  sensitivity  of  the  two-dimensional  runoff  model  to  the  spatial  variability  of  precipitation 
intensity  is  explored.  Both  stochastically  generated  and  polarimetric  weather  radar  estimated  static 
rainfall  fields  are  input  to  the  runoff  model  for  differing  rainfall  durations  and  spatial  resolutions. 
The  sensitivity  of  the  runoff  model  to  rainfall  intoisity  spatial  variability  is  obse^ed  to  decrease 
with  increasing  rainfall  duraticm  for  all  spatial  resoluticms  tested.  The  effect  of  precipitation  data 
resolution  is  observed  by  comparistm  with  results  from  the  finest  resolution  input  data.  The 
statistical  deviations  in  model  perftnmance  due  to  resolution  coarsening  are  presented. 

The  effect  of  storm  motion  on  the  outflow  calculated  the  two-dimensional  runoff  nxxlel  is 

observed  using  rectangular,  constant  intensity  precipitation  fields.  The  precipitation  fields  are 
moved  in  16  directions  over  the  watershed,  at  ten  speeds  per  direction.  A  dimensionless  storm 
speed  and  stcnm  direction  is  identified  which  prodiKes  the  largest  effect  on  the  magnitude  of  the 
outflow  hydrograph  peak.  This  dimensionless  storm  velocity  is  on  die  order  of  the  runoff  length 
divided  by  the  time  to  equiUbrium  for  an  equivalent  stationary  storm. 

Real-time  polarimetric  weather  radar  estimated  rainfall  fields  are  input  to  the  two-dimensional 
runoff  model  in  simulations  with  and  without  infiltration  on  one  32  km^  and  one  121  knfi 
watershed  to  determine  the  effect  of  input  data  spatial  resolution  on  model  performance.  In 
simulations  without  infiltraticm,  the  sn^er  watershed  is  more  sensitive  to  input  data  resolution 
because  of  increased  variability  in  rainfall  volume.  Statistical  deviations  fix>m  finest  resolution 
results  are  presented.  In  simulatitms  widi  infiltration,  the  interaction  between  the  scale  of  the  input 
rainfall  data  and  infiltration  processes  in  ctxisidetable.  As  input  rainfall  data  resolution  is 
coarsened,  excess  rainfall  amounts  decrease  substantially,  indicating  that  two-dimensional, 
distributed  parameter,  physically  based  mnoff  models  must  be  calilxated  for  one  input  data 
resoluticMi.  If  the  resolution  of  Ae  precipitaticm  data  is  changed,  the  model  must  be  re-calibrated. 
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INFERENCE  OF  HORIZONTAL  TEMPERATURE  GRADIENTS 
USING  PASSIVE  RADIOMETRIC  METHODS 

Peter  Olsson 

Atmospheric  Science  Department 


There  exist  many  situations  in  nature  whoe  relatively  strong  horizontal  temperature  gradients 
are  present  in  the  boundary  layer.  The  purpose  of  this  work  is  to  investigate  the  impact  of 
horizontal  temperature  gradients  on  the  in&aied  radiance  properties  of  the  boundary  layer  and 
devise  a  scheme  for  inferring  the  gradient  magnitude  from  these  radiance  properties. 

The  temperature  and  spectral  dependence  of  the  radiance  of  various  atmospheric  constituents 
are  examined  and  the  IR  portion  of  the  spectrum  SOOcm-i  ^  v  ^  SOOcm-^  is  shown  to  be  most 
sensitive  to  temperature  induced  radiance  changes.  The  two  impcHtant  radiating  gases  at  these 
wavenumbers  are  CO2  and  H2O.  The  spectral  radiative  properties  of  these  constituents  are 
discussed.  The  variability  of  the  water  vapor  mixing  ratio  is  shown  to  be  an  important  factOT  in 
this  remote  sensing  application.  A  model  is  discuss^  which  numerically  estimates  the  radiance  of 
the  boundary  layer  in  Ae  presence  of  a  hcxizcHital  tenq)erature  gradient  The  results  of  this  model 
demonstrate  the  possibili^  of  estimating  the  gradient  magnitude  from  narrow  band  azimuthally 
scanned  radiance  measurements. 

Two  parameters,  the  attenuation  length  and  the  centered  nramalized  radiance,  are  introduced 
and  their  relationship  to  the  gradient  magnitude  is  explored.  Using  these  parameters,  a  method  is 
developed  which  permits  the  inference  of  temperature  gradient  magnitude  from  infrared  radiance 
measurements  and  local  pressure,  temperature  and  mixing  ratio  values.  The  success  of  this 
technique  is  demonstrated  by  the  accurate  recovery  of  gradient  magnitudes  frxxn  calculated  radiance 
data.  Finally,  consideration  is  given  to  sources  of  error  and  uncertainty  in  die  measurement 
process  and  the  inqiact  of  these  on  the  inference  of  the  gradient  magnitude. 
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THE  EFFECT  OF  MOVING  RAINSTORMS  ON  OVERLAND  FLOW 
USING  ONE-DIMENSIONAL  FINITE  ELEMENTS 

Jerry  R.  Richardson 
Civil  Engineering  Department 


Rainfall  events  are  rarely  stationary.  That  is  to  say  that  over  any  watershed,  stonns  are 
spatially  and  tempraally  variable.  Unfortunately,  most  rainfall-runoff  noodels  usually  ctxisider 
precipitation  to  be  stationary  and  uniftnm.  However,  evidence  in  the  literature  indicates  that  the 
movement  of  storms  can  have  a  significant  influence  on  runoff  hydrographs.  Therefore,  a  need  for 
numerical  models  which  can  accommodate  the  movement  of  precipitation  is  recognized.  This 
dissertation  provides  a  better  understanding  of  the  influence  of  block  moving  precipitation  on 
surface  runoff  hydrographs.  Two  one-dimensional  finite  element  models  were  developed  and  used 
for  these  investigations. 

Analytical  expressions  for  each  partial  derivative  of  the  momentum  equaticm  are  developed  fw 
stationary  stonns.  For  moving  rain  storms,  the  behavior  of  die  terms  of  the  momentum  equation  is 
sho^  to  be  very  complex.  However,  bas^  on  the  analysis  of  the  momentum  equatitxi  for 
stationary  and  block  moving  steams,  it  is  found  that  only  the  fitU-dynamic  and  kinematic  forms  of 
the  momentum  equation  are  suitable  for  simulation  of  overland  flow.  The  kinematic  approximation 
is  suitable  for  ovo’land  flow  simulatitxi  for  most  cases  of  stationary  and  moving  stonns.  The  time 
to  peak  of  moving  storms  was  shown  to  depend  on  the  traverse  time  of  storms  and  can  be  several 
times  greater  than  the  time  to  equilibrium  for  stationary  storms.  Fot  partial-equilibrium 
hydrographs,  the  maximum  peak  dischai:ge  of  equivalent  block  moving  storms,  was  shown  to 
occur  when  the  dimensionless  storm  velocity  was  equal  to  0.3  far  laminar  flow  and  0.5  for 
turbulent  flow  using  Maniting's  equation.  The  maximum  peak  discharge  occurs  when  the  flood 
wave  and  the  stt>im  front  arrive  at  the  downstream  end  of  Ae  plane  at  the  same  time.  Using  the 
results  of  these  studies,  and  the  finite  element  model,  success^  simulation  of  block  moving 
storms  over  a  laboratory  watershed  was  conducted.  These  simulations  indicate  that  the  one¬ 
dimensional  finite  element  model  can  be  used  to  simulate  runoff  from  nooving  stonns  over  a  system 
of  planes  and  channels. 
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THE  BROADBAND  RADIATIVE  PROPERTIES  OF  CIRRUS  CLOUDS 
DEDUCED  FROM  AIRCRAFT  MEASUREMENTS  DURING  FIRE 

William  L.  Smith,  Jr, 

Atmospheric  Scirace  Department 


The  bulk  radiative  and  microphysical  pr(^)erties  of  five  cirms  clouds  sampled  via  the  NCAR 
Sabreliner  on  four  days  during  the  FIRE  ^t  cirrus  BFO  are  described.  These  cirrus  systems, 
which  developed  under  a  variety  of  sync^dc  weather  condidons,  occurred  at  various  aldtudes  and 
ranged  in  geometric  thickness  from  about  2.0  to  4.5  km.  A  bro^band,  infrared  radiadve  transfer 
model  is  employed  to  deduce  the  intact  of  the  cirrus  layers  m  infrared  radiadon.  This  ttKxlel 
isolates  the  effect  of  the  atmospheric  gases  fix>m  that  of  the  cloud  ice  water  pomitdng  retrieval  of 
the  cloud  emittance  (Ccu)  and  profiles  of  the  mass  absorpdon  coefficient  (iO* 

Fot  the  five  cirrus  cloud  cases,  the  total  cloud  emittance,  €cu,  ranged  frtxn  about  0.4  to  0.8  and 
the  deduced  emittance  profiles  appear  as  similar  funcdons  of  ice  water  path  (IWP).  Furthermore, 
the  mass  absorpdon  coefficient,  K,  is  found  to  decrease  with  increasing  pardcle  size  ranging  from 
about  0.48  in  the  top  of  one  layer  to  about  0.007  m^g-i  near  the  base  of  anmher.  This 
reladonship  is  somewhat  dssimilar  from  one  cirrus  system  to  the  next  suggesting  the  significant 
effect  of  some  uiuneasured  microphysical  property.  Small  particles,  which  have  been  shown  by 
other  authors  to  be  prevalent  in  cirrus  clouds  via  die  spectral  characterisdcs  renootely  sensed  in  the 

8-12  M-m  window  region,  are  a  likely  suspect  Broadband,  infrared  absorpdon  coefficients  (a)  are 
also  computed  and  found  to  exhibit  a  sin^ar  temperature  dependence  as  data  recendy  presented  by 
other  authors. 

The  horizontal  variabilides  in  the  shwtwave  and  infrared  properdes  of  these  cirrus  systems  are 
explored.  The  range  of  variadon  in  the  shortwave  prc^ierties  are  found  to  be  similar  to  the 

observed  range  in  the  infrared.  Good  correladon  was  found  between  the  shortwave  albedo  (p)  and 
upward  effecdve  emittance  (€*  T ).  A  scatter  plot  of  these  two  parameters  agreed  well  with 
theoredcal  calculadons  assuming  an  asymmetry  parameter  of  0.7.  Downward  effecdve  emittances 
(€*  i)  were  foimd  to  range  from  about  0.4  to  0.8,  while  the  shortwave  effecdve  exdncdon  (Q 
ranged  from  0  to  0.45.  e*  i  and  ^  were  not  well  correlated  owing  to  cloud  heterogeneides. 

Finally,  the  current  state  of  cirrus  radiadon  parameterizadons  was  briefly  assessed  in  reladon  to 
this  d^  set  and  there  appears  to  be  sufficient  observadonal  evidence  to  support  the  inidal 
development  of  parameterizadon  schemes  for  general  circuladon  and  climate  models. 
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MECHANICS  OF  HEADCUT  MIGRATION  IN  RILLS 

Otto  R.  Stein 

Civil  Engineering  Department 


The  proposed  conceptual  model  relates  two-dimensional  headcut  migration  to  sediment 
detachment  just  upstream  and  just  downstream  from  a  headcut  If  upstream  erosion  dominates,  the 
headcut  tends  to  obliterate  itsefr  as  it  migrates  upstream,  eventually  looming  indistinguishable 
from  the  eroding  channel's  bed  slope.  If  downstream  erosion  dc»ninates,  the  headcut  face  erodes 
from  below  and  a  definable  headcut  with  a  near  vertical  face  migrates  upstream  with  time. 

A  criterion  to  determine  which  migration  mode  occurs  is  formulated  using  both  dimensional 
analysis  and  by  equating  hydraulic  and  sediment  detachment  equations.  Dimensional  analysis 
results  in  a  time  s^e  ratio  of  upstream  to  downstream  erosion  related  to  the  upstream  and 
downstream  sediment  detachment,  bed  slope,  Reynolds  number  and  the  drq)  number,  which 
represents  the  dimensionless  headcut  drop  height  A  physically  based  analysis  of  hydraulics  and 
SMliment  detachment  yields  a  relationship  between  the  dimensionless  terms.  The  resulting  equaticm 
for  headcut  stability  is  favorably  compart  with  a  total  of  eleven  laboratory  measurements  of 
headcut  migration  on  cohesive  soil. 

The  maximum  scour  depth  and  the  total  volume  of  eroded  material  produced  by  the  impinging 
jet  just  downstream  from  a  headcut  are  analyzed.  A  previously  developed  method  for  the 
pre^ction  of  the  ultimate,  or  equilibrium,  scour  depth  in  non-cdliesive  bed  material  is  modified  ro 
include  scour  from  cohesive  bed  material.  This  method  equates  the  sediment  detachment  potential 
of  the  bed  material  to  the  diffusion  of  a  jet 

The  same  approach  is  used  to  determine  the  time  rate  of  the  maximum  scour  depth. 
EHmensional  analysis  relates  the  ratio,  maximum  scour  depth  at  any  time  over  the  predicted  ultimate 
scour  depth,  to  the  jet  Reynolds  number,  Froude  number  and  a  dimensionless  time  scaled  to  jet 
properties.  The  change  in  scour  depth  is  analytically  shown  to  proceed  at  two  distina  rates.  For 
some  period  of  time  Som  the  initiation  of  scour,  scour  rate  is  independent  of  time  because  the  bed 
is  within  the  jet  potential  core  and  diffusion  has  not  reduced  the  maximum  jet  veloci^.  Beyond 
this  time  period  jet  diffusion  decreases  the  rate  of  scour  and  this  rate  decrees  with  increasing  time 
and  scour  depth.  In  the  limit  the  predicted  ultimate  scour  depth  is  approached.  The  dimensional 
and  physically  based  analyses  for  scour  depth  ctxnpare  favorably  with  experimental  data.  This 
data  includes  measurements  of  the  ultimate  scour  (tepth  and  the  change  in  scour  depth  with  time  for 
ten  runs  on  cohesive  soil,  eight  runs  on  sand  dso  =  1.5mm  and  six  runs  on  sand  dso  =  0.15mm. 
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MICROWAVE  REMOTE  SENSING  OF  A  PRECIPITATING  ATMOSPHERE 


Joseph  Turk 

Electrical  Engineering  Department 


Throughout  the  10-100  GHz  spectrum,  the  scattering  and  absorption  properties  of  atmospheric 
hydrometeors  and  gases  vary  widely.  As  a  result,  multifi%quency  techniques  for  passive 
microwave  estimation  of  precipitation  have  been  proposed.  The  emission-based,  lower  frequency 
observations  (<  37  GHz)  effectively  respond  at  lower  altitudes,  while  the  scattering-based  higher 
frequencies  ^  37  GHz)  respond  at  higher  altitudes,  thus  providing  a  potential  means  to  link  self- 
consistent  precipitation  estimates  to  the  multifrequency  tightness  temperatures  Tb.  Unlike  many 
schemes  which  link  Tb  to  the  ’iquid  water  content,  this  study  emphasizes  the  use  of  scattering- 
biised  channels  to  infer  iiiformation  on  the  ice  water  content  The  high  altitude  ice  region  remains 
fairly  unobscured  to  a  spacebome  sensor,  suggesting  the  implementation  of  the  scattering-based 
channels  in  a  top-down  methodology  as  the  firet  step  in  a  retrieval  of  the  cloud  vertical  stracture. 

Extensive  model  simulations  of  the  upwelling  Tb  were  performed  throughout  clouds  whose 
bulk  microphysical  properties  were  deduced  from  multiparameter  radar  data  and  cloud  model 
output  obtained  during  the  1986  Cooperative  Huntsville  Meteorological  Experiment  (COHMEX). 
A  highly  accurate  plane-parallel  microwave  radiative  transfer  model  with  Nfre  phase  matrices  was 
used  to  solve  the  equation  of  transfer  in  a  scattering  atmosphere  with  a  land  surface.  COHMEX 
aircr^  radiometer  measurements  at  18, 37,  and  92  GHz  were  compared  with  the  output  of  the 
multiparameter  radar-initialized  radiative  transfer  model.  Deviations  at  92  GHz  were  attributed  to 
uncertainties  in  the  bulk  density  and  size  distribution  of  the  ice  region.  Using  both  multiparameter 

radar  observations  and  cloud  model  output,  the  37-85  GHz  Tb  difference  (ATb)  was  found  to  be 
sensitive  to  the  amount  of  integrated  ice  water  path  (IWP)  lying  above  the  rain,  independent  of  ice 

density.  Presence  of  coexisting  cloud  water  and  ice  had  the  net  effect  of  compressing  the  ATb  for 
a  given  IWP.  Over  convective  regions,  85  GHz  Tb  depressions  <  200  K  were  noted  to  be 
relatively  insensitive  to  the  amount  of  underlying  liquid  water  content,  and  mainly  dependent  upon 
the  optical  thickness  of  the  ice  layer. 
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SOCIAL  LOAFING  AND  VIGILANCE  DECREMENT  WHILE  MONITORING 
DIFFERENT  LEVELS  OF  SIMULATED  AIR  TRAFFIC 

Paul  A.  Weiler 
Psychology  Department 


Using  a  computerized  simulation  of  an  air  traffic  control  display,  this  study  examined  the 
phenomenon  of  social  loafing  (a  decrease  of  individual  contribution  in  a  group  setting)  and 
vigilance  decrement  (the  decrease  in  monitoring  performance  over  time).  A  total  of  96  introductory 
psychology  students  served  as  subjects  ina2X2X2X3  split-plot  design.  The  independent 
variables  were  the  difficulty  of  the  task  (easy,  hard),  perceived  uniqueness  of  an  individual's 
contribution  (unique,  non-unique),  and  i^ntifiability  of  individual  performance  (identifiable,  non- 
identifiable).  Additionally,  three  lO-minute  time  blocks  within  the  simulation  were  examined 
Groups  of  four  subjects  were  instmcted  on,  and  then  participated  in,  a  3()-minute  simulation  of  an 
air  traffic  controller  radar  screen.  Subjects  were  asked  m  monitor  the  simulated  traffic  for  a  critical 
stmtilus  of  '999'  in  the  display.  Subject  performance  was  recorded  as  the  latency  of  response  to  a 
critical  stimulus  once  present^  Results  showed  significant  main  effects  for  task  difficulty  and  for 
drw  block  within  the  simulation.  Additional  results  showed  the  failure  of  the  identifiability  and 
uniqueness  manipulations  hypothesized  as  necessary  for  social  loafing  to  occur.  Both  failures  may 
be  related  to  the  lack  of  perception  of  a  group  setting  by  the  subjects.  It  is  recommended  that  future 
studies  increase  subject  perception  of  group  membership. 
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THE  INFLUENCE  OF  TOPOGRAPHY  ON  COLORADO 
FRONT  RANGE  SNOWSTORMS 

Douglas  A.  Wesley 
Atmospheric  Science  Department 


This  study  utilizes  both  an  extensive  set  of  observations  and  mesoscale  model  simulations  to 
isolate  and  describe  the  important  influences  of  complex  terrain  on  Colorado  Front  Range  winter 
storms,  with  an  emphasis  on  snowfall  distributions.  Specifically,  the  interaction  of  various  types 
of  cold,  low-level  air  masses  with  topography  and  the  larger-sc^e  flow  is  described.  Frequently, 
the  heaviest  snowfall  does  not  coincide  with  the  steepest  terrain  gradients,  as  might  be  expected, 
due  to  this  interaction.  Field  measurements  of  several  snowstorms  during  the  1987-89  time  period 
include  special  CLASS  (Cross-chain  Loran  Atmospheric  Sounding  System)  vertical  profiles, 
standard  National  Weather  Service  (NWS)  and  National  Meteorological  Center  (NMC)  siuface  and 
rawinsonde  da^  surface  observations  collected  by  a  trained  snow-spotter  network,  profiler  data, 
Doppler  reflectivity  and  velocity  scans,  and  Geostationary  Operation^  Environmental  Satellite 
(GOES)  images.  The  numerictti  predictions  are  produced  by  the  CSU  Regional  Atnoospheric 
Modeling  System  (RAMS),  using  two-  and  three-dimensional  non-hydrostatic  simulations.  The 
versions  of  tiie  mo^l  employed  in  this  study  utilized  both  horizontally  homogeneous  initializations 
and  initial  fields  made  up  of  NMC  upper-air  grids  and  rawinsondes.  Most  simulations  en^loyed 
full  microphysics.  A  omit  comparison  betv'een  nxxlel-predicted  dynamic  and  microphysical  fields 
and  observational  data  is  described  for  the  30-31  March  1988  snowstorm  along  the  Colorado  Front 
Range. 

Of  particular  interest  in  this  study  is  the  role  of  trapped  cold  air  masses  over  the  foothills  and 
adjacent  plains  during  the  evolution  of  snow-prodwting  synoptic  scale  disturbances.  Specifically, 
these  include  two  types  of  {vecipitating  easterly  ("upslq)e")  flows  along  the  Colorado  Front 
Range;  cold-air  daring  situations  and  arctic  outlmalu.  By  no  means  are  the  dynamical  features 
of  these  two  types  of  storms  mutually  exclusive.  The  investigation  includes  detailed  case  studies 
of  each  storm  type.  The  cold  air  daintning  process  ccmcentrates  snowfall,  typically  in  the  foothills 
and  adjacent  areas  to  the  east  A  blocked  low-level  stable  layer  causes  overrunning  in  a  north- 
south  band  over  and  near  the  foothills.  A  quite  different  scenario,  the  cold  outbre^  requires  the 
establishment  of  an  arctic  air  mass  in  the  low  levels.  In  this  case,  overrunning  occurs  as  moist 
westerly  winds  aloft  flow  over  the  arctic  air  mass,  producing  snowfall  in  the  foothills  and  plains. 

The  problems  with  forecasting  precipitation  amounts  during  these  situaticms  are  well-known. 
This  project  isolates  two  of  the  primary  mechanisms  which  have  caused  many  of  these  problems. 
Both  of  these  dynanucal  processes  must  be  considered  by  the  wintertime  forecaster  as  major 
contributors  to  heavy  snowfall  in  this  region. 
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3.  CONCLUSION 


The  Geosciences  Fellowship  program  has  been  very  successful  firan  the  perspective  of  the 
Fellows  and  from  the  research  projects  that  conqjrise  the  Center  for  Geosciences.  Nineteen 
graduate  degrees  had  been  granted  by  Colcaado  State  University  (CSU)  as  a  result  of  the  program 
up  to  the  time  that  this  Final  Repeat  was  written.  It  is  expected  that  the  Fellowship  program  will 
have  contributed  at  least  partially  to  26  degrees  by  the  time  that  Fellows  who  still  are  working  on 
their  programs  of  study  have  finished. 

The  Fellowship  program  also  has  been  very  successful  in  maintaining  and  strengthening  the 
Geosciences  research  community  through  the  en^loyment  of  its  graduates.  The  following 
graduates  are  en^)loyed  by  the  Center  fw  Geosciences  as  professicmal  scientists: 

Robert  McCoy 
Joseph  Turk 

Others  are  employed  at  fetteral  research  laboramries,  including: 

James  Bossert  <  Los  Alamos  Nsdonal  Lab,  Department  of  Energy 
Raul  Alvarez  -  Environmental  Protection  Agency 

Private  industry  has  employed  several  of  the  former  Fellows: 

Thaddeus  Cline 
Russell  Huonder 
Jeny  Richardson 
Wilikm  Smith 
Paul  Weiler 

Former  Fellow  Otto  Stein  has  become  an  Assistant  ProfessOT  at  Montana  State  University.  The 
largest  group  of  Fellows  are  working  on  additional  advanced  degrees  at  Colorado  State  University 
and  elsewhere: 


Deborah  Anthony  -  CSU 
Charlotte  Atwater  -  CSU 
Keely  Costigan  -  CSU 
Andrew  Jones  -  CSU 
Michael  Meyers  -  CSU 
Glen  Moglen  -  MTT 
Peter  Olsson  -  CSU 

The  Fellows  contributed  significantly  to  the  research  activities  of  the  Center  for  Geosciences. 
The  interdisciplinary  character  of  research  at  the  Center  was  partcularly  enhanced  by  the  efforts  of 
Uie  Fellows.  For  example,  Joseph  Turk  worked  extensively  with  Andrew  Jones  on  the  theory  and 
int^retation  of  attrx)^hetic  microwave  measurements  and  with  Fred  Ogden  on  using  radar  ^ta 
as  input  to  a  hydrologic  model.  The  Fellows  also  were  senior  authors  or  co-authors  on  a  number 
of  scientific  publications. 


